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We  investigate  theoretically  and  experimentally  the  statistical  properties  of  dc  current  through  an  open 
quantum  dot  subject  to  ac  excitation  of  a  shape-defining  gate.  The  symmetries  of  rectification  current  and 
photovoltaic  current  with  respect  to  applied  magnetic  field  are  examined.  Theory  and  experiment  are  found  to 
be  in  good  agreement  throughout  a  broad  range  of  frequency  and  ac  power,  ranging  from  adiabatic  to  nona- 
diabatic  regimes. 

DOI:  10.1103/PhysRevB.71.241309  PACS  number(s):  73.23. -b,  73.40.Ei,  73.50.Pz 


Transport  in  mesoscopic  systems  subject  to  time-varying 
fields  combines  elements  of  nonequilibrium  physics  and 
quantum  chaos.  This  combination  extends  the  scope  of  me¬ 
soscopic  physics  and  is  likely  to  be  important  in  quantum 
information  processing,  where  fast  gating  and  quantum  co¬ 
herence  are  both  required.  Of  particular  importance  is  the 
ability  to  control  external  fields  applied  to  the  mesoscopic 
system  and  to  distinguish  effects  of  these  fields  on  quantum 
dynamics  of  the  system.  For  example,  two  distinct  contribu¬ 
tions  to  direct  current  through  an  open  quantum  dot  due  to  an 
oscillating  perturbation  have  been  identified1'2  and  observed 
experimentally  in  Ref.  3. 

In  this  paper,  we  investigate  the  statistical  properties  of  dc 
currents  resulting  from  an  applied  ac  electric  field  over  a 
wide  range  of  excitation  frequencies,  paying  particular  atten¬ 
tion  to  the  presence  or  absence  of  symmetry  with  respect  to 
magnetic  field  in  various  regimes.  Theoretical  analysis  is 
based  on  recently  developed  time-dependent  random  matrix 
theory.4'5  Experiments  use  a  gate-defined  GaAs  quantum  dot 
subject  to  ac  excitation  of  a  gate  at  MHz  to  GHz  frequencies. 
At  low  excitation  frequencies  rd]  (rd  is  the  electron 
dwell  time  in  the  dot)  the  present  theoretical  results  are  con¬ 
sistent  with  those  obtained  by  adiabatic  approximations.1,6 
However,  the  analysis  is  applicable  over  a  wider  range  of 
frequencies  hu>  ^Et,  where  ET=hl  rcross  is  the  Thouless  en¬ 
ergy  and  Tcross  is  the  electron  crossing  time  of  the  dot.  At 
hut's. ET,  the  system  may  be  studied  by  methods  developed 
for  disordered  bulk  conductors.7-8 

Three  distinct  contributions  to  dc  current  through  the  dot 
can  be  identified,  resulting  from;  (i)  an  applied  dc  bias;  (ii) 
an  ac  bias  at  the  excitation  frequency  (i.e.,  rectification 
effects2);  and  (iii)  photovoltaic  effects.5,9  Although  the  ori¬ 
gins  of  the  rectification  and  photovoltaic  effects  are  different, 
these  two  effects  are  hard  to  distinguish  in  experiment  with 
rectification  current  being  dominant.10  The  purpose  of  the 
present  Rapid  Communication  is  to  describe  distinct  features 
between  the  rectification  and  photovoltaic  effects.  We  restrict 
our  attention  to  a  one-parameter  excitation,  noting  that  while 
in  the  adiabatic  regime  one-  and  two-parameter  excitations 
affect  the  system  differently,1  beyond  the  adiabatic  regime, 
to  Sr)1,  the  differences  disappear.5 

The  Hamiltonian  of  electrons  in  the  dot  in  the  presence  of 
a  magnetic  flux  $  is  represented  by  a  Hermitian  MXM  ma¬ 
trix  //ft)  =  7T,|,+V  cos  tot.  We  assume  that  electron  dynamics 


in  the  dot  is  fully  chaotic.11  Then  the  time-independent  part 
Htj,  may  be  considered  as  a  random  realization  of  a.  MXM 
matrix  from  a  Gaussian  unitary  ensemble  with  the  mean 
level  spacing  S{  and  M~ET/SV  Perturbation  V  is  a  matrix 
from  a  Gaussian  orthogonal  ensemble  characterized  by  the 
strength  C0=  it  TtV1 1 M1 .  The  parameter  C0  determines  the 
energy  displacement  of  an  electron  state  due  to  the  applied 
perturbation  V.  The  contact  between  the  left  (right)  lead  and 
the  dot  contains  N\(Nr)  open  channels,  we  enumerate  chan¬ 
nels,  a,  in  the  left  {a=  1 , . ..  ,Nj)  and  the  right  (a=N\ 
+ 1 , ...  ,Ach)  contacts,  Ncb=Ni+Nr.  The  corresponding  ex¬ 
perimental  setup  is  shown  in  Fig.  1. 

The  dc  current  7'1’  through  the  dot  is  determined  by  the 
scattering  matrix  [5 <b{t,t')\ap  (see  Ref.  5) 
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is  the  distribution  function  of  electrons  in  channel  a  at  tem¬ 
perature  T  and  voltage  Vjt).  At  sufficiently  low  frequencies 
u><iEclh  (Ec  is  the  dot  charging  energy)  VJ t)  is  simply 
related  to  the  bias  V(t)  across  the  dot:  Va{t)  =  AaaV(t).  Ele¬ 
ments  of  the  diagonal  matrix  A  are  A aa=Nr/Nch  for  l^a 
^Nh  and  A aa=-NfNch  for  N,<a^Nch. 

We  consider  the  bias  V(f)  across  the  dot  in  the  form 
y(t)  =  V0+Vulcos(a>t+  tp^).  The  dc  current  through  the  dot  to 
first  order  in  dc  bias  V0  and  ac  bias  V(t>  is12 

P  =  7*+/f  +  ^Vo,  7f  =  fX,  (3) 

where  the  first  term  represents  the  photovoltaic  current  T6 
=P(Va=0),  see  Eqs.  (1)  and  (2).  The  second  and  third  terms 
in  Eq.  (3)  represent  the  contributions  to  the  current  due  to  dc 
bias  VQ  and  ac  bias  VM,  respectively. 
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FIG. _ L  (Color  online)  Symmetry  factor  Sk 

=  {SG\ ®  SGk'1’)  /  ((SGf)2)  as  a  function  of  frequency  w  for  k= 0  (up¬ 
per  panel)  and  k=  1  (lower  panel)  at  two  values  of  temperature  T 
and  power  Co  of  the  ac  excitation.  Inset:  micrograph  of  the  device 
and  a  schematic  picture  of  applied  voltages. 
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and  =  0  (see  Ref.  6).  In  Eqs.  (6)  and  (7)  the  angle 

brackets  (•  ■  ■)  stand  for  the  averaging  with  respect  to  realiza¬ 
tions  of  T-iq,.  Functions  and  B^,.t  describe  the  distribu¬ 
tion  function15  of  electrons  in  the  dot  in  the  presence  of 
time-dependent  electric  fields  and  kernels  K±=K~t,.gj  de¬ 
scribe  the  evolution  of  electron  states16  in  these  fields.  Both 
functions  BU'\  and  B]f\  and  the  kernels  K~,:  „  contain  the 
diffuson  IX 1 1  ,f2,  r)  =  exp[-/Jir(r, t)r/t]  or  the  Cooperon 
C(r1,T2,f)  =  exp[-(y)/JT(r,f)c/r].  Here,  F  (T,t)  =  ye+ yv 
+4 C0  sin2wf  sin2(wr/2),  where  ye=  SlNcil/2Tr  is  the  electron 
escape  rate  and  yv  is  the  electron  phase  relaxation  rate  due  to 
inelastic  processes. 

In  the  experiment,  the  ac  bias  V M  results  from  capacitive 
coupling  between  the  leads  and  the  gate  on  which  the  ac 
voltage  is  applied  (see  the  inset  in  Fig.  1).  Therefore  VM  is 
proportional  to  the  amplitude  of  the  ac  voltage  at  the  gate. 
Assuming  that  Cn  is  linear  in  the  applied  power  to  the  gates, 
we  write  V(J  =  «fiJyC(l/  ye.  The  coefficient  aM  has  units  of  volt¬ 
age  and  is  independent  of  realizations  of  the  quantum  dot 
(we  disregard  fluctuations  of  V  over  different  realizations  of 
Hcj,).  Therefore,  the  correlation  functions  of  the  rectification 
current  Tf =gfVw  are  determined  by  the  correlators  of  SGf 

(if’Tf)  =  SGf).  (8) 


where  Go=NiNr/ Nc^  is  the  classical  conductance  and 
8Gk(t)  =  Jq^" SGk(t) codt / 277  stands  for  time  averaging  of  the 
“instantaneous  conductance”  (k= 0,  1) 

SGf(t)  =  {  dTFk(r)  J  ddcos(k(a>6  +  cpk)) 

XTrjA^tf-^A^tf+^J, 

7 ikBTTHi  2irkBT  sin(coT/2) 

°  sinh(7 TkgTrUi)  1  hco  sinh(  irk BTr/h) 

(5) 

We  observe  that  SGl(t)  =  SGQ(t)cos(u)t  + ipf)  in  the  adiabatic 
limit  hu>< max{Ach<5i ,kBT}  considered  in  Refs.  2,13. 

Below  we  study  the  variance  of  the  photovoltaic  current 
Jfh  and  the  conductances  and  g'f  with  respect  to  random 
realizations  of  the  Hamiltonian  X,[,.  Following  Refs.  5  and 
14,  we  find  in  the  limit  A'ch  =s>  1  and  at  magnetic  fields  <f> 
§><l>0\Wch/M  destroying  the  weak  localization  ( <F0=hc!e ) 


We  also  notice  that  in  the  limit  A'ch  §>  1  the  correlation  func¬ 
tion  of  the  photovoltaic  current  /A  and  the  rectification  cur¬ 
rent  7f  vanishes:6  (7^  7f)  =  0. 

First  we  use  Eqs.  (7)  and  (8)  to  analyze  the  magnetic  field 
symmetry  of  the  rectification  current  if.  Although  in  the 
adiabatic  limit  a><yjh  the  rectification  current  ff  is  sym¬ 
metric  with  respect  to  magnetic  field  inversion  (<& ->-<!>),  at 

higher  frequencies  co^yjfi  the  symmetry  of  If  is  sup¬ 
pressed.  Indeed,  the  magnetic  field  symmetry  is  related  to  the 
time  inversion  symmetry.  For  a  harmonic  field  at  frequency 
oj,  the  time-inversion  symmetry  holds  only  on  time  scales 
much  smaller  than  1  /«.  Transport  through  the  system  is  de¬ 
termined  by  times  of  the  order  of  hi  ye  and  consequently  the 
magnetic  field  symmetry  of  the  rectification  current  breaks  if 
kaye.  We  plot  the  ratio  =  SG~®) I  ((SGf)2)  as  a 

function  of  hu>!  ye  in  Fig.  1.  .S’ j  =  1  represents  the  symmetric 
current  If  «  SGf  with  respect  to  magnetic  field  inversion.  In 
the  adiabatic  regime  this  symmetry  originates  from  the  On- 
sager  symmetry17  of  the  dc  conductivity  (see  Eq.  (5)  and 
Refs.  2,3).  As  the  frequency  increases,  .S’ ,  vanishes,  signal- 
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ling  the  suppression  of  the  magnetic  field  symmetry.  There¬ 
fore,  the  absence  of  magnetic  field  symmetry  no  longer 
serves  as  a  feature  which  allows  one  to  distinguish  the  pho¬ 
tovoltaic  current  7^h  from  the  rectification  current  if. 

We  notice  that  the  magnetic  field  symmetry  of  the  dc 
conductance  <5G|f  is  more  sturdy  than  the  symmetry  of  the 
rectification  current  (see  Fig.  1).  Particularly,  at  temperatures 
kBT a  ye,  dc  conductance  cSG®  is  nearly  symmetric  at  fre¬ 
quencies  ha>^kBT,  since  the  dc  correlation  function  is  deter¬ 
mined  by  processes  on  a  time  scale  h/kBT.  The  symmetry  is 
not  fully  suppressed  even  at  w>kBT/h\  the  suppression  de¬ 
pends  on  C0/ye. 

We  apply  Eqs.  (6)-(8)  to  the  analysis  of  the  experiment.3 
The  quantum  dot  used  in  the  experiment  has  an  area  A 
~0.7  /iiii2.  Relevant  energy  scales  are  the  Thouless  energy 
Et~  160  /ueV  and  the  mean  level  spacing  bj  =  ttjc 
— 10  /xeV.  The  measurements  were  performed  at  the  base 
electron  temperature  T— 200  mK  (kBT=5Aye).  From  the 
size  of  conductance  fluctuations  without  ac  fluctuation  on  the 
gate,  we  estimate  the  dephasing  rate  y^—0.2 ye  (see  Ref.  18 
for  details)  and  thus  disregard  yv  and  energy  relaxation  rate 
y,,  y  in  our  quantitative  analysis. 

In  the  upper  panel  of  Fig.  2  we  show  the  variance  of  the 
conductance  as  a  function  of  the  incident  power  P  at  W2tt 
=  5.56  GHz  (ha>l  ye~l  .2).  We  also  plot  the  variance  of  the 
conductance  calculated  from  Eq.  (7)  (k= 0)  at  temperature 
T=5AyJ  kB.  Assuming  that  the  ratio  C0/ye  is  proportional  to 
the  power  P  of  the  ac  excitation  applied  to  the  gate,  i.e., 
C0/ye=P/P0,  we  rescale  P  to  obtain  the  best  fit  of  the  ex¬ 
perimental  points  by  the  curve  of  Eq.  (7).  We  find  P0= 9 
X  1(T8  W. 

In  the  lower  panel  of  Fig.  2  we  show  the  correlators 
(7+'1,7±"’)  of  the  measured  current.  Although  the  traces  of  the 
magnetic  field  sweeps  (see  Fig.  3  in  Ref.  3)  look  quite  asym¬ 
metric  for  the  measured  current,  the  antisymmetric  correlator 
(l+'l,J-"r)  is  not  significantly  smaller  than  its  symmetric  coun¬ 
terpart.  We  notice  however,  that  if  the  averaging  is  per¬ 
formed  over  n  realizations,  the  measured  correlator  (7+*7_4>) 
can  be  estimated  as  (/+<l7+<1')/  xn  (n  —  50  in  the  experiment). 

We  plot  the  variance  of  the  photovoltaic  current,  Eq.  (6), 
as  a  function  of  Q/ ye  for  kBT=  5.4ye  and  hu>=l ,2ye  (to 
facilitate  numerics,  we  used  the  approximation  ha>9>ye).  We 
emphasize  that  the  horizontal  shift  between  the  data  points 
and  curve  is  fixed  by  the  fit  in  the  upper  panel  for  the  dc 
conductance  and  there  are  no  fitting  parameters  for  the  vari¬ 
ance  of  the  current  (along  the  vertical  axis).  At  C0S  ye,  the 
variance  of  the  measured  current  changes  quadratically  in 
C0/  ye,  consistent  with  quadratic  dependence  on  C0  of  the 
theoretical  curves  for  (flh)  and  therefore  our  assumption  that 
C0  is  proportional  to  the  power  of  the  ac  excitation  is  justi¬ 
fied.  At  C()£  yc  the  variance  of  the  measured  current  starts 
saturating.  This  saturation  is  expected  for  large  power  as¬ 
ymptote  of  the  photovoltaic  current  due  to  the  spreading  of 
the  distribution  function  of  electrons  in  the  dot.5  Some  de¬ 
viation  between  the  experimental  points  and  the  theoretical 
curve  is  expected  due  to  the  approximation  /Vch  S>  1  used  for 
derivation  of  Eq.  (6)  (Nch=  2  in  the  experiment). 

For  illustration,  we  also  plot  the  correlation  functions  of 
the  rectification  current,  using  Eq.  (8)  for  <p ,  =  0  and  a(il 
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FIG.  2.  (Color  online)  Upper  panel:  Variance  of  the  conductance 
as  a  function  of  the  ac  excitation  power  P=P0C0I  ye  at  (oIItt 
=  5.56  GHz  (hco= 7.2ye)  and  the  theoretical  result  of  Eq.  (7)  with 
k= 0.  We  use  PQ= 9X  10-8  W  and  kBT=5Aye.  Lower  panel:  Sym¬ 
metric  (solid  triangles)  and  antisymmetric  (open  triangles)  current 
correlators  as  a  function  of  ac  excitation  strength  C0.  The  solid  line 
shows  the  variance  of  the  photovoltaic  current  Eq.  (6)  with  param¬ 
eters  fixed  by  the  fit  in  the  upper  panel.  The  dashed  and  dotted  lines 
show  the  symmetric  and  antisymmetric  correlators  of  the  rectifica¬ 
tion  current  Eq.  (8)  with  aot=0A5fi.a>/e  and  <Pi  =  0. 

=  0A5hu>le.  For  the  rectification  current,  the  saturation  at 
large  power  is  not  expected:  according  to  Fig.  2, 

(C0/  ye)a  with  a  —  0.6. 

We  similarly  discuss  the  data  for  u>/2tt=2A  GHz.  Per¬ 
forming  the  fit  of  the  experimental  values  of  the  conductance 
fluctuations  and  the  result  of  Eq.  (7)  with  k= 0  and  tempera¬ 
ture  T=5Aye/ kB,  we  find  the  relation  between  the  strength  of 
the  perturbation  Q  and  the  power  P=P0Cq/ ye  with  P0 
=  2.5  X  10-7  W. 

In  Fig.  3  we  show  the  symmetric  and  antisymmetric  cur¬ 
rent  correlators  for  2.4  GHz.  For  comparison  we  plot  by  a 
solid  line  the  variance  of  the  photovoltaic  current  7jfh,  calcu¬ 
lated  from  Eq.  (6)  at  ha>=  3.1  ye.  We  observe  that  the  fluctua¬ 
tions  of  the  measured  current  significantly  exceed  (by  a  fac¬ 
tor  —  100)  the  expected  magnitude  for  the  photovoltaic 
current,  and  therefore  are  likely  due  to  the  rectification  of  the 
bias  across  the  dot.  The  low  power  data  can  be  fitted  by  Eq. 
(8)  with  aw=A.lhu>i e. 

The  above  choice  for  a01'z{fi u>,kBT}/e  limits  the  applica¬ 
bility  of  the  linear  expansion  Eq.  (3)  to  small  powers  of  the 
ac  excitation  C0,  such  that  C0/ ye^{ha>/ eaj2.  The  higher 
order  corrections  in  the  bias  Va  do  not  restore  magnetic  field 
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FIG.  3.  (Color  online)  Symmetric  (solid  triangles)  and  antisym¬ 
metric  (open  triangles)  current  correlators  at  to!2rn=2A  GHz  as  a 
function  of  power  P  =  P0Cq/ ye,  with  P0=2.5X  10-7  W.  The  solid 
line  shows  the  variance  of  the  photovoltaic  current  Eq.  (6)  at  tem¬ 
perature  T=5Aye/ks  and  frequency  u>=3.1  ye/fi.  The  dashed  and 
dotted  lines  show  the  symmetric  and  antisymmetric  correlators  of 
the  rectification  current  Eq.  (8)  with  am=4.7ft.col e  and  <pi=0. 

symmetry,  which  is  in  apparent  contradiction  to  the  observed 
symmetry  of  the  measured  current  at  larger  powers  (at 
Co/ ye~  1  in  Fig.  3).  We  attribute  the  restoration  of  magnetic 
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field  symmetry  to  dephasing  due  to  dot  heating  by  the  dissi¬ 
pative  current.  Increasing  the  power  P  at  fixed  oj  drives  the 
system  into  the  adiabatic  regime  since  the  heating  makes  the 
ratio  ha>/ (ye+  yv)  decrease.  As  shown  in  Fig.  1,  the  rectifi¬ 
cation  current  is  symmetric  in  the  adiabatic  regime.  The  as¬ 
sumption  that  yv  increases  as  power  P  increases  is  consistent 
with  the  observed  change  of  the  correlation  field  for  the  cur¬ 
rent  fluctuations.  Indeed,  according  to  Fig.  4  in  Ref.  3,  as  the 
power  changes  from  10  to  104nW,  ye  +  y^  increases  by  fac¬ 
tor  of  4  and  becomes  larger  than  hoo.  In  this  regime  the 
rectification  current  is  mostly  symmetric  with  respect  to 
magnetic  field  inversion  (see  Fig.  1). 

In  summary,  we  studied  ensemble  fluctuations  of  dc  cur¬ 
rent  through  an  open  quantum  dot  subject  to  oscillating  per¬ 
turbation.  We  showed  that  as  frequency  of  the  perturbation 
increases,  magnetic  field  symmetry  of  the  current  disappears, 
regardless  of  the  mechanism  of  the  current  generation.  We 
demonstrated  that  the  power  behavior  of  the  current  fluctua¬ 
tions  is  an  important  tool  to  distinguish  effects  of  an  ac  ex¬ 
citation  on  dc  current. 
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